Formation of unusual Cr5+ charge state in CaCr0.5Fe0.5O3 perovskite
Dai Jian-Hong1, 2, Zhao Qing1, Sun Qian1, Zhang Shuo3, Wang Xiao1, 2, Shen Xu-Dong1, 2, Liu Zhe-Hong1, 2, Shen Xi1, 2, Yu Ri-Cheng1, Chan Ting-Shan4, Li Lun-Xiong5, Zhou Guang-Hui6, Yang Yi-feng1, 2, †, Jin Chang-Qing1, 2, Long You-Wen1, 2, ‡
Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
School of Physics, University of Chinese Academy of Sciences, Beijing 100049, China
Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201203, China
“National” Synchrotron Radiation Research Center, 101 Hsin-Ann Road, Hsinchu 30076, Taiwan, China
Institute for Brain Research and Rehabilitation, South China Normal University, Guangzhou 510631, China
Department of Physics and Synergetic Innovation Center for Quantum Effects and Applications of Hunan, Hunan Normal University, Changsha 410081, China

 

† Corresponding author. E-mail: yifeng@iphy.ac.cn ywlong@iphy.ac.cn

Abstract
Abstract

A new oxide CaCr0.5Fe0.5O3 was prepared under high pressure and temperature conditions. It crystallizes in a B-site disordered Pbnm perovskite structure. The charge combination is determined to be Cr5+/Fe3+ with the presence of unusual Cr5+ state in octahedral coordination, although Cr4+ and Fe4+ occur in the related perovskites CaCrO3 and CaFeO3. The randomly distributed Cr5+ and Fe3+ spins lead to short-range ferromagnetic coupling, whereas an antiferromagnetic phase transition takes place near 50 K due to the Fe3+–O–Fe3+ interaction. In spite of the B-site Cr5+/Fe3+ disorder, the compound exhibits electrical insulating behavior. First-principles calculations further demonstrate the formation of charge combination, and the electron correlation effect of Fe3+ plays an important role for the insulting ground state. CaCr0.5Fe0.5O3 provides the first Cr5+ perovskite system with octahedral coordination, opening a new avenue to explore novel transition-metal oxides with exotic charge states.

1. Introduction

Transition-metal oxides with unusual valence states show fascinating physical properties. For example, an exotic charge-ordering transition that involves competing dimer and trimer orderings is observed in the recently discovered compound Fe4O5 obtained by high-pressure synthesis methods.[1,2] Both temperature and pressure induced intersite charge transfers are found to occur in the A-site ordered quadruple perovskite LaCu3Fe4O12 with the presence of unusually high Cu3+ and Fe3.75+ valence states, accompanying with a series of sharp variations in unit cell, magnetism, and electrical transport properties as well as negative thermal expansion.[35] Charge disproportionation of Fe4+ takes place in the distorted perovskite CaFeO3 , leading to metal-to-insulator and orthorhombic-to-monoclinic structural phase transitions.[6,7] On the other hand, in the cubic perovskite SrFeO3, the high Fe4+ state is stable and there is not any charge variation at different temperatures. However, the antiferromagnetic (AFM) ordering of SrFeO3 is very sensitive to chemical doping. A small amount of Fe4+ substituted by another high Co4+ or Ni4+ can lead to strong ferromagnetic (FM) interactions with high Curie temperature above room temperature.[810]

As for the element of chromium, the valence states such as Cr3+ with half-filled orbitals and Cr6+ with empty 3d orbitals exist widely. Although the Cr4+ state ( can occur in the simple binary oxide CrO2 under ambient pressure, high pressure is necessary to prepare ternary ACr4+O3 perovskites such as CaCrO3, SrCrO3, and PbCrO3, where interesting physical properties are observed. The crossover from local to itinerant electron features gives rise to anomalous electronic behaviors for the orthorhombic CaCrO3 and cubic SrCrO3 perovskites.[1115] Large volume collapse and Pb–Cr intermetallic charge transfer emerge in PbCrO3.[1618] By comparison, the chromium oxide with Cr5+ is rare to know. Though the zircon-type RCr5+O4 ( earth and Y) compounds with tetrahedral coordinated CrO4 units can form by adopting low-temperature wet chemical methods, the oxygen is readily to release on heating.[1923] At higher temperatures, these compounds always decompose to Cr3+-based RCrO3 perovskites. At present, there is no report for the discovery of Cr5+-involved perovskite oxide with octahedral coordinated Cr5+O6 unit ( .

In this paper, we report a new perovskite CaCr0.5Fe0.5O3 (CCFO) synthesized at 8 GPa and 1473 K. Since the well-studied CaCrO3 and CaFeO3 respectively possess quadrivalent Cr4+ and Fe4+ at room temperature,[6,7,11,12,24,25] one may think, at first glance, that the Cr and Fe ions in CCFO should also possess similar charge states (i.e., ). In sharp contrast, however, our experimental results indicate the presence of Cr5+ and Fe3+ states. CaCr0.5Fe0.5O3 thus provides the first example with unusually high Cr5+ state in octahedral coordinated environment of perovskite.

2. Experimental and calculation details

The solid-state reaction synthesis of polycrystalline CCFO was carried out using a cubic-anvil-type high pressure apparatus. Highly pure ( ) CaO, Fe2O3, CrO2, and CrO3 powders were used as the starting materials. These reactants with a stoichiometric mole ratio of 4:1:1:1 were mixed thoroughly in an agate mortar within an argon-filled glovebox. The mixed powders were then sealed into a platinum capsule with 2.8 mm diameter and 4.0 mm length. The capsule was treated at 8 GPa and 1473 K for 30 min in the high-pressure apparatus. When the heating process was finished, the sample was quenched to room temperature and then the pressure was slowly released to ambient pressure. The crystal quality and structure of CCFO were examined by powder x-ray diffraction (XRD) using a Huber diffractometer with Cu radiation at 40 kV and 30 mA. The XRD data were collected in the angle (2θ) range from 10° to 100° with steps of 0.005° and analyzed by the Rietveld refinement using the GSAS program.[26] Thermogravimetric (TG) measurement was performed using a Setaram TG-DTA system in mixed H2 (4%) and N2 (96%) gas flow on heating from 300 K to 1473 K with a speed of 5 K/min. The selected area electron diffraction (SAED) was carried out using a Philips-CM200 field emission transmission electron microscopy operated at 200 kV. The Fe and Cr K- and L-edge x-ray absorption spectroscopy (XAS) data were collected at beamline 14W of the Shanghai Synchrotron Radiation Facility and beamline 16A1 of the “National” Synchrotron Radiation Research Center of Taiwan. The dc magnetic susceptibility (χdc) and magnetization were measured using a Quantum Design superconducting quantum interference device magnetometer. The ac magnetization ( ) and resistivity were measured on a physical property measurement system (Quantum Design, PPMS-9T).

Density functional theory calculations were performed using the full-potential linearized augmented plane-wave method implemented in WIEN2K code with the generalized-gradient approximation Perdew-Burke-Ernzerhof (GGA-PBE) exchange-correlation energy.[27] The muffin-tin radii were set to 2.1 a.u. for Ca, 1.70 a.u. for O, and 1.90 a.u. for Cr/Fe. The maximum modulus for the reciprocal vectors was chosen such that and we used 1000 k-point meshes for the Brillouin zone.

3. Results and discussion

Figure 1 shows the XRD pattern of the high-pressure product CCFO. The Rietveld refinement shows that CCFO crystallizes into an ABO3-type perovskite structure with orthorhombic space group Pbnm, where the B-site Cr and Fe are randomly distributed at the special atomic site 4b(0, 0.5, 0), forming oxygen six-fold coordinated Cr/FeO6 octahedra. To further confirm the disordered Cr/Fe distribution, high-resolution SAED was performed along different characteristic axes like [1-10] zone axis (not shown here). We did not find any trace for long-range B-site ordering in CCFO. This observation is essentially different from the theoretical prediction on the analogous oxide Sr2CrFeO6, which is claimed to be a B-site ordered double perovskite composed of high-spin Cr4+ and low-spin Fe4+ charge states.[28] Table 1 lists the refined structural parameters as well as the satisfied goodness factors for the Rietveld analysis of CCFO.

Fig. 1. (color online) X-ray diffraction pattern and structure refinement results of CaCr0.5Fe0.5O3. The observed (circles), calculated (line), and difference (bottom line) are shown. The top ticks indicate the allowed Bragg reflections of CCFO with space group Pbnm. The bottom ticks indicate the impurity phase of CaFe2O4. The two-phase refinement shows that the amount of impurity phase is 4.1 wt.%.
Table 1.

Refined structure parameters of CaCr0.5Fe0.5O3 at room temperature. Space group: Pbnm; atomic positions: Ca (x, y, 0.25), Cr/Fe (0, 0.5, 0), O1 (x, y, , O2 (x, y, 0.25). Occup.: occupancy factor.

.

To confirm the oxygen content in CCFO, the temperature dependence of TG measurement was carried out. As shown in Fig. 2, the compound is very sensitive to heating and starts to lose some weight above 350 K, which may reflect the metastable state of the Cr5+/Fe3+ charge combination with the presence of unusually high Cr5+ valene state as described later. According to the TG loss between 350 K and 1473 K as well as the residual products (CaO, Fe, and CaCr2O4), the oxygen content is determined to be 2.97 ± 0.02. This value is very close to the stoichiometric oxygen composition for CCFO.

Fig. 2. Temperature dependence of thermogravimetric loss for CaCr0.5Fe0.5O3.

We now turn to the determination of the charge states of Fe and Cr in CCFO by hard XAS K-edge measurements. Figure 3(a) shows the Fe K-edges of CCFO together with those of Fe2O3 and SrFeO3, which are used as standard Fe3+ and Fe4+ references with oxygen octahedral coordination, respectively. The chemical shift of the Fe K-edge defined at –0.8 of the normalized intensity can be used to identify the Fe valence state.[29,30] We find that the absorption edge of CCFO shifts to lower energy by 1.5 eV relative to that of SrFeO3, but at the same energy as that of Fe2O3, revealing the Fe3+ valence state of CCFO. The Cr K-edge of CCFO is shown in Fig. 3(b), where Cr2O3 and SrCrO3 are used as Cr3+ and Cr4+ references with octahedral coordination, respectively. Unfortunately, at this stage there is no octahedral coordinated Cr5+ reference that can be applied for comparison. Anyway, in Fig. 3(b) one can see that the Cr K-edge shifts to higher energy by 1.3 eV from Cr3+ oxide Cr2O3 to Cr4+ oxide SrCrO3 and by another 1.0 eV from SrCrO3 to CCFO. Taking into account the confirmed Fe3+ state as well as the nearly stoichiometric chemical composition of CCFO, the systematical energy shift observed in the Cr K-edge strongly suggests the formation of an unusual Cr5+ state fulfilling the charge balance requirement for the charge combination.

Fig. 3. (color online) The hard XAS of CaCr0.5Fe0.5O3 measured at room temperature for (a) the Fe-K edge, and (b) the Cr-K edge. The spectra of some related references are also shown for comparison. The rectangles show the normalized intensity –0.8, which can be used to judge the valence state.

The unusual Cr5+/Fe3+ charge combination in CCFO can be further confirmed by soft XAS at the Fe- and Cr-L2,3 edges. As shown in Fig. 44(a), the energy position of the Fe-L2,3 edges in CCFO is remarkably lower than that of SrFe4+O3 by about 0.8 eV, ruling out the formation of Fe4+. By comparison, the Fe-L2,3 XAS spectrum of CCFO situates at similar energy and has comparable spectral feature with that of LaFeO3, confirming the presence of high-spin Fe3+ valence state in CCFO. Since there is no Cr5+ reference with an octahedral coordination, the Cr-L2,3 XAS of Cr2O3 and SrCrO3 were measured as Cr3+ and Cr4+ references with octahedral coordination, respectively. As shown in Fig. 4(b), with increasing valence state from Cr3+ in Cr2O3 to Cr4+ in SrCrO3, the main peak positions of both L3 and L2 edges shift toward higher energy by ∼1.0 eV. As the absorption spectrum of CCFO is concerned, it also experiences a similarly systematic energy shift, but displays substantially different spectral features from that of the Cr4+ reference SrCrO3 probably due to the enhanced p–d hybridization between Cr and the coordinating oxygen atoms, suggesting that the unusually high Cr5+ state instead of Cr4+ occurs in CCFO, in agreement with the K-edge measurements mentioned above. Both the hard and soft XAS results thus demonstrate the charge combination. The Cr5+/Fe3+ charge states presented in the current CCFO greatly differ with the Cr4+ state observed in CaCrO3 and Fe4+ in CaFeO3. The newly synthesized oxide CCFO provides a rare example possessing an unusual Cr5+ state in a perovskite structure with Cr5+O6 octahedral coordination.

Fig. 4. (color online) The soft XAS of CaCr0.5Fe0.5O3 measured at room temperature for (a) the Fe-L2,3 edges, and (b) the Cr-L2,3 edges. The spectra of some related references are also shown for comparison.

The unexpected Cr5+/Fe3+ charge states of CCFO are also consistent with the magnetism as characterized by both dc and ac magnetization measurements. Figure 5(a) shows the temperature dependence of χdc measured at 0.1 T. With decreasing temperature to , the χdc sharply increases, whereas there is no visible anomaly in ac magnetization (Fig. 5(b)). Moreover, the zero-field-cooling (ZFC) and field-cooling (FC) χdc curves start to separate from each other below TSR. These behaviors may imply the presence of some short-range FM interactions below TSR due to the random distribution of Fe3+ and Cr5+ spins. On further cooling to , both the ZFC and FC χdc curves experience a remarkable kink, suggesting an AFM phase transition. When the ac magnetization is measured at different frequencies (see Fig. 5(b)), we find that the AFM cusp occurring around 50 K is not dependent on the measurement frequency at all, confirming that the AFM ordering should be moderately long-range in nature. Above 200 K, the inverse χdc can be well fitted by the Curie–Weiss law (Fig. 5(a)), producing the Curie constant and the Weiss temperature . The negative Weiss temperature is consistent with the low-temperature long-range AFM ordering, and the considerably large value of θ compared with TN may imply some spin frustrations caused by the competing FM and AFM interactions mentioned above. Accordeing to the Curie constant, the effective magnetic moment of CCFO is calculated to be /f.u., which is roughly comparable to the spin-only theoretical value for charge conbination (4.36 /f.u).

Fig. 5. (color online) Temperature dependence of (a) dc magnetic susceptibility and its inverse measured at 0.1 T, and (b) ac magnetization ( ) measured at different frequencies for CaCr0.5Fe0.5O3. (c) Field dependence of magnetization measured at different temperatures. The purple line in (a) shows the Curie–Weiss fitting.

Figure 5(c) shows the dc magnetization measured at different temperatures. Above TSR, the linear magnetization behavior is coherent with the paramagnetism. Below TSR, however, one can find a small amount of magnetic hysteresis due to the short-range FM coupling as well as the possible canted antiferromagnetism (below TN). In CCFO, there exist two distinct magnetic ions (i.e., Cr5+ and Fe3+). According to the Goodenough-Kanamori rules,[31,32] the Cr5+–O–Fe3+ and Fe3+–O–Fe3+ interaction pathways with the average band angle ∼159° (see Table 1) will generate FM and AFM superexchange interactions, respectively. Our first-principles calculations also show the FM interaction for the 3d1 electrons of Cr5+ via the Cr5+–O–Cr5+ superexchange pathways (not shown here). Consequently, both Cr5+–O–Fe3+ and Cr5+–O–Cr5+ interactions are responsible for the short-range FM correlation of CCFO, whereas the AFM ordering should arise from the Fe3+–O–Fe3+ superexchange interaction.

Because of the strong electronic correlated effects, the 3d5-electron perovskite systems of Fe3+ such as RFeO3 usually exhibit Mott insulating behaviors.[33] For the 3d1-electron systems like LaTiO3 and YTiO3 with distorted perovskite structure, the electrical insulating behavior is also observed,[3436] whereas the undistorted cubic 3d1 perovskite SrVO3 displays metallic conductivity.[37] In the current CCFO perovskite oxide, there exists considerable structural distortion. The compound also shows the insulating behavior, as featured by the increasing resistivity on cooling as well as the well fitted 3D Mott variable-range-hopping model (see Fig. 6 and the inset), although the B-site Fe3+ (3d5) and Cr5+ (3d1) are distributed disorderly. Moreover, as seen from Table 1, the Cr-O distances in a single Cr/FeO6 octahedron in CCFO can be differed by about 4.5%. This significant difference may suggest possible orbital ordering for the 3d1 electrons of Cr5+ ions, as observed in another 3d1 system Sr2VO4 with tetragonal distorted perovskite structure.[38]

Fig. 6. (color online) Resistivity as a function of temperature for CaCr0.5Fe0.5O3. The red line in the inset shows the fitting result on the basis of the 3D Mott variable-range-hopping model between 200 K and 300 K.

Since the detailed electronic properties of Cr5+ in perovskite structure with octahedral coordination are completely absent to date, it is very interesting to further investigate through first-principles calculations based on the density functional theory (DFT). Considering that the DFT is unfavorable to treat a completely disordered system, for simplify, we assume some ordered crystal structures to gain insights into the properties of the experimental CCFO. The electronic structures of CCFO were preliminary calculated for different hypothetical Cr/Fe orderings such as the G-type, A-type, and C-type. By comparison, the G-type (i.e., rocksalt type) Cr/Fe ordered structure is optimized in system energy and we obtain the lattice parameters a = 5.123 Å, b = 5.375 Å, and c = 7.580 Å with space group P21/n. In this structure model, we obtain an A-type AFM spin ground state, with a zero total magnetic moment. The band structures and the partial density of states (DOS) of the magnetic ground state were calculated with Ueff = 3 eV for Cr and 4 eV for Fe, following the common choices in the literature,[3941] and the results are shown in Fig. 7, indicating an insulating behavior with an energy gap of about 0.3 eV. Both are consistent with the experimental observation in CCFO. Our calculated moment is about 1.0 μB/Cr and 4.1 μB/Fe, so there is only one d-electron on each Cr ion and a fully polarized and half-filled d-shell on each Fe ion. We therefore conclude that the Cr ions have a valence of +5, while the Fe ions have a valence of +3 in this hypothetical G-type Cr/Fe ordered perovskite structure. In sharp contrast, for other types of Cr/Fe orderings mentioned above, we obtain only Cr4+ ions. In addition, in the mentioned G-type Cr/Fe ordered structure, if we do not consider the Coulomb interaction for Cr5+ and Fe3+ both, a metallic ground state instead of an insulating one as observed in experiment is obtained. However, one can still obtain an insulating ground state (albeit with a much reduced gap) if Uef = 0 is set only for Cr5+, suggesting that the Coulomb interaction on Fe3+ ions is crucial for the insulating conductivity of CCFO.

Fig. 7. (color online) First-principles numerical results for the band structures and the partial density of states of the hypothetical G-type Cr/Fe ordered double perovskite Ca2CrFeO6.
4. Conclusion

In summary, by using high-pressure and high-temperature synthesis conditions, we obtained a new perovskite oxide CaCr0.5Fe0.5O3 with nearly stoichiometric oxygen content. Although previous theoretical calculation predicted the analogous compound Sr2CrFeO6 to be a B-site ordered double perovskite with Cr4+/Fe4+ charge combination, our results reveal that CCFO crystallizes in a simple ABO3 perovskite structure with randomly distributed Cr and Fe ions at the B site. Fe- and Cr-K and L edge XAS measurements suggest that the charge states are Fe3+ and Cr5+ in CCFO, providing the first example that shows an unusual Cr5+ state in perovskite structure with Cr5 +O6 octahedral coordination. In accordance with the Cr5 +/Fe3+ disordered perovskite structure, the random Cr5 +–O–Fe3+ and Cr5+–O–Cr5+ interactions cause some short-range FM couplings below ∼120 K. On the other hand, the Fe3+–O–Fe3+ superexchange interaction generates an AFM ordering at a lower temperature 50 K. The resistivity data of CCFO follows the Mott 3D variable-range-hopping mechanism, indicating the insulating conductivity. Our first-principle theoretical calculations based on the hypothetical G-type Cr/Fe ordered structure further illustrate the formation of Cr5+/Fe3+ charge combination. The electorn correlation effect of Fe3+ is closely related to the insulating ground state of CCFO.

Reference
[1] Ovsyannikov S V Bykov M Bykova E Kozlenko D P Tsirlin A A Karkin A E Shchennikov V V Kichanov S E Gou H Y Abakumov A M Egoavil R Verbeeck J McCammon C Dyadkin V Chernyshov D van Smaalen S Dubrovinsky L S 2016 Nat. Chem. 8 501
[2] Lavina B Dera P Kim E Meng Y Downs R T Weck P F Sutton S R Zhao Y 2011 Proc. Natl. Acad. Sci. USA 108 17281
[3] Long Y W Hayashi N Saito T Azuma M Muranaka S Shimakawa Y 2009 Nature 458 07816
[4] Long Y W Kawakami T W Chen T Saito T Watanuki T Nakakura Y Liu Q Q Jin C Q Shimakawa Y 2012 Chem. Mater. 24 2235
[5] Long Y W Shimakawa Y 2010 New J. Phys. 12 063029
[6] Takano M Nakanishi N Takeda Y Naka S Takada T 1977 Mat. Res. Bull. 12 923
[7] Woodward P M Cox D E Moshopoulou E Sleight A W Morimoto S 2000 Phys. Rev. B 62 844
[8] Kawasaki S Takano M Takeda Y 1996 J. Solid State Chem. 121 174
[9] Long Y W Kaneko Y Ishiwata S Tokunaga Y Matsuda T Wadati H Tanaka Y Shin S Tokura Y Taguchi Y 2012 Phys. Rev. B 86 064436
[10] Seki H Hosaka Y Saito T Mizumaki M Shimakawa Y 2016 Angew. Chem. Int. Ed. 55 1360
[11] Zhou J S Jin C Q Long Y W Yang L X Goodenough J B 2006 Phys. Rev. Lett. 96 046408
[12] Long Y W Yang L Lv Y Liu Q Jin C Zhou J Goodenough J B 2011 J. Phys.: Condens. Matter 23 355601
[13] Komarek C Möller T Isobe M Drees Y Ulbrich H Azuma M Fernández-Díaz M T Senyshyn A Hoelzel M André G Ueda Y Grüninger M Braden M 2011 Phys. Rev. B 84 125114
[14] Ortega-San-Martin L Williams A J Rodgers J Attfield J P Heymann G Huppertz H 2007 Phys. Rev. Lett. 99 255701
[15] Lee K W Pickett W E 2009 Phys. Rev. B 80 125133
[16] Arévalo-Lópezá M Dossantos-García A J Alario-Franco 2009 Inorg. Chem. 48 5434
[17] Cheng J Kweon K E Larregola S A Ding Y Shirako Y Marshall L G Li Z Y Li X dos Santos A M Suchomel M R Matsubayashi K Uwatoko Y Hwang G S Goodenough J B Zhou J S 2015 Proc. Natl. Acad. Sci. USA 112 1670
[18] Yu R Hojo H Watanuki T Mizumaki M Mizokawa T Oka K Kim H Machida A Sakaki K Nakamura Y Agui A Mori D Inaguma Y Schlipf M Rushchanskii K Z Ležaić M Matsuda M Ma J Calder S Isobe M Ikuhara Y Azuma M 2015 J. Am. Chem. Soc. 137 12719
[19] Jiménez E Isasi J Sáez-Puche R J 2000 Alloy. Compd. 312 53
[20] Tsirlin A A Rabie M G Efimenko A Hu Z Saez-Puche R Tjeng L H 2014 Phys. Rev. B 90 085106
[21] Sáez-Puche R Jiménez E Isasi J Fernández-Díaz M T García-Muñoz J L 2003 J. Solid State Chem. 171 161
[22] Long Y W Liu Q Lv Y Yu R Jin C 2011 Phys. Rev. B 83 024416
[23] Long Y W Yang L X Yu Y Li F Y Yu R C Jin C Q 2007 Phys. Rev. B 75 104402
[24] Bhobe P A Chainani A Taguchi M Eguchi R Matsunami M Ohtsuki T Ishizaka K Okawa M Oura M Senba Y Ohashi H Isobe M Ueda Y Shin S 2011 Phys. Rev. B 83 165132
[25] Takano M Nasu S Abe T Yamamoto K Endo S Takeda Y Goodenough J B 1991 Phys. Rev. Lett. 67 3267
[26] Larson A C and Von Dreele R B Report No. LAUR 86-748 (Los Alamos National Laboratory Los Alamos NM 1994)
[27] Blaha P Schwarz K Madsen G K H Kvasnicka D Luitz J 2013 WIEN2K: An Augmented Plane Wave + Local Orbitals Program for Calculating Crystal Properties Karlheinz Schwarz, Technische Universitat Wien Wien, Austria 2013
[28] Lee K W Ahn K H 2012 Phys. Rev. B 85 224404
[29] Croft M Sills D Greenblatt M Lee C Cheong S W Ramanujachary K V Tran D 1997 Phys. Rev. B 55 8726
[30] Poltavets V V Croft M Greenblatt M 2006 Phys. Rev. B 74 125103
[31] Goodenough J B 1955 Phys. Rev. 100 564
[32] Kanamori J 1959 J. Phys. Chem. Solids 10 87
[33] Xu W M Naaman O Rozenberg G K Pasternak M P Taylor R D 2001 Phys. Rev. B 64 094411
[34] Pavarini E Biermann S Poteryaev A Lichtenstein A I Georges A Andersen O K 2004 Phys. Rev. Lett. 92 176403
[35] Ulrich C Gössling A Grüninger M Guennou M Roth H Cwik M Lorenz T Khaliullin G Keimer B 2006 Phys. Rev. Lett. 97 157401
[36] Khaliullin G Maekawa S 2000 Phys. Rev. Lett. 85 3950
[37] Lan Y C Chen X L He M 2003 J. Alloys Compd. 354 95
[38] Eremin M V Deisenhofer J Eremina R M Teyssier J van der Marel D Loidl A 2011 Phys. Rev. B 84 212407
[39] Streltsov S V Korotin M A Anisimov V I Khomskii D I 2008 Phys. Rev. B 78 054425
[40] Feng H L Arai M Matsushita Y Tsujimoto Y Guo Y Sathish C I Wang X Yuan Y H Tanaka M Yamaura K 2014 J. Am. Chem. Soc. 136 3326
[41] Deng H S Liu M Dai J H Hu Z W Kuo C Yin Y Y Yang J Y Wang X Zhao Q Xu Y J Fu Z M Cai J W Guo H Z Jin K J Pi T Soo Y Zhou G H Cheng J G Chen K Ohresser P Yang Y F Jin C Q Tjeng L H Long Y W 2016 Phys. Rev. B 94 024414